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Long, cylindrical metal nanowires have recently been observed to form and be stable for seconds 
at a time at room temperature. Their stability and structural dynamics is well described by a 
continuum model, the nanoscale free-electron model, which predicts cylinders in certain intervals of 
radius to be linearly unstable. In this paper, I study how a small, localized perturbation of such 
an unstable wire grows exponentially and propagates along the wire with a well-defined front. The 
front is found to be pulled, and forms a coherent pattern behind it. It is well described by a linear 
marginal stability analysis of front propagation into an unstable state. In some cases, nonlinearities 
of the wire dynamics are found to trigger an invasive mode that pushes the front. Experimental 
procedures that could lead to the observation of this phenomenon are suggested. 

PACS numbers: 47.54.-r, 66.30.Pa, 68.65.La, 



Front propagation into unstable states occurs in many 
areas of physics, chemistry, and biology (see Ref. [l| for a 
recent review), and is often related to pattern formation 
mechanisms. In this paper, I show that metal nanowires, 
whose dynamics can be described by a continuum model, 
the nanoscale free-electron model (NFEM) 0, [1, Q , ex- 
hibit such front propagation into unstable cylinders. 

Recent transmission electron microscopy experiments 
[S H 0, HI have observed gold and silver nanowires to 
form long cylinders, with diameters of order one nanome- 
ter, that are stable for seconds at a time at room temper- 
ature. A theoretical description within the NFEM shows 
[3, Q that such self-assembly is natural if conditions are 
such that the atom mobility allows the wire to explore 
its configuration space and reach its equilibrium shape, 
as is the case at room temperature for these metals. 

Due to the high number of surface atoms — with low 
coordination numbers — in such nanowires, surface ef- 
fects are particularly important, and favor wire break- 
up due to the Rayleigh instability [l(| EH ■ It has been 
shown, using the NFEM [J, E3 , that the quantum con- 
finement of electrons in the cross-section of the wire pro- 
vide electron-shell effects — similar to those well-known in 
cluster physics [l3| — that compete with surface effects, 
and stabilize cylindrical wires for a finite range around 
magic radii [ll|, [12], as well as a number of number of 
wires with broken axial symmetry [l4l . Il5j . 

The NFEM @, H, [| is a continuum model where the 
atomic structure is replaced by a uniform, positively 
charged background, and the emphasis is put on the elec- 
tronic structure. An extension of the model [3j includes 
ionic dynamics through surface self-diffusion, which is ex- 
pected to dominate the structural dynamics in such thin 
wires. One of its important predictions is that a random 
wire naturally evolves into a universal equilibrium shape 
consisting of a perfect cylinder of a magic radius, con- 
nected to thicker leads 3, 9]. Recently, the rich kink dy- 
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namics of the NFEM has been described and shown 
to be qualitatively similar to the thinnin g m echanism of 
gold nanowires observed experimentally [17| . proving it 
to be a suitable, if simplified, model of the structural 
dynamics of metal nanowires. 

This paper is concerned with radii outside of the inter- 
vals of stability, and studies their dynamics under surface 
self-diffusion. A localized perturbation to an unstable 
cylinder is found to grow exponentially, and to propa- 
gate along the unstable wire with a well defined front, 
which is a coherent pattern forming front, and can be 
either pulled or pushed depending on the wire radius. 

The paper is organized as follows: In Sec.|H the NFEM 
is introduced and some of its main results directly rele- 
vant to the present article are summarized. Numerical 
simulations of the dynamics of an unstable cylinder are 
presented in Sec. UH while the front propagation anal- 
ysis is developed in Sec. IIIII Section IIVI discusses the 
results and some experimental setups that could detect 
front propagation into unstable metal nanowires. 



I. THE NANOSCALE FREE-ELECTRON 
MODEL 

The NFEM is a continuum model of open metallic 
nanosystems with an emphasis on the electronic struc- 
ture, which is treated exactly 0,3. It is thus particularly 
suitable as a model of metal nanowires, and successfully 
describes many of their equilibrium 0, E2, E2, EH, EM E3| 
and dynamical [!, 0, H, EE H3| properties in simple phys- 
ical terms. 

The discrete atomic structure is replaced by a uniform, 
positively-charged background (Jellium), which provides 
a confining potential for electrons. Electronic degrees of 
freedom are described using a free-electron model, thus 
neglecting interactions, except inasmuch as they rescale 
macroscopic quantities such as the bulk energy density 
ujb, and the surface tension a s 0, EH EH- This model is 
particularly suitable for simple metals, with good screen- 
ing and a close-to-spherical Fermi surface, such as those 
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with a single s-electron conduction band at the Fermi 
surface. Such conditions are fulfilled for alkali metals 
like sodium, and for noble metals such as gold and silver, 
although d electrons may play a role for noble metals. 

While the NFEM has no such restriction, only axisym- 
metric wires are considered in this paper. This choice 
greatly simplifies the numerical treatment of the dynam- 
ics, as the wire shape can be described by a single radius 
function R{z,t). It is justified by the facts that (i) the 
most stable wires are axisymmetric [lj], HH , and (ii) the 
dynamics tends to decrease surface area, and thus further 
favors axisymmetric wires. 

A nanowire being an open system, the electronic en- 
ergy is given by the grand canonical potential fl e . Like 
any extensive thermodynamic quantity, f2 e can be writ- 
ten as a Weyl expansion [2l| — a series in geometrical 
quantities such as system volume V and surface area S — 
complemented by a mesoscopic, fluctuating contribution 
Sfl: 



n e [R(z)} = u B v + a s s + sn, 



(1) 



where the values of lob and as niay be chosen to match 
the bulk properties of the metal to be described. As 
results are independent of ujb, its free-electron value 
ujb = —2EFkp/15ir 2 is used, while the surface tension 
is set to a s = 1.256 N/m, a value appropriate for the 
description of gold (22j. 

Assuming the wire cross section varies slowly along the 
wire (adiabatic approximation), the mesoscopic contribu- 
tion may be written as 



5Sl[R(z)] = / dzV sheU [R{z)}, 
Jo 



(2) 



where the electron-shell potential V s heii(R) can be com- 
puted using a semi-classical approximation [4] . V s heii (R) , 
depicted in Fig. [T] (top panel) for a cylindrical wire as a 
function of its radius R, is responsible for stabilizing wires 
of magic radii, which correspond to its deep minima. 

The ionic dynamics is taken to be classical and can be 
assumed to occur mainly through surface self-diffusion, as 
most atoms in thin metal wires are surface atoms 0, [23j . 
The evolution equation for the radius function R(z, t) 
derives from ionic mass conservation 



a 

V a ot oz 



(3) 



where V a = 'in 2 /k^ F is the volume of an atom, and the 
z-component J z of the total surface current is given by 
Fick's law: 



J, = - 



p s D s 2wR(z,t) dp 
k B T ^/T+Jd z Mj 2 d^' 



(4) 



Here ps and D$ are, respectively, the surface density 
of atoms and the surface self-diffusion coefficient, and 
d z R = dR/dz. 



The precise value of D$ is not known for most met- 
als, but it can be removed from the evolution equation 
by rescaling time to the dimensionless variable r = u>ot, 
with the characteristic temperature-dependent frequency 
= PsDsTp/T. For comparison to experimental time 
scales, one can estimate that for quasi-one-dimensional 
diffusion D s rs vua? exp(— Eg/ksT), where vu is the De- 
bye frequency, a is the lattice spacing, and E s is an acti- 
vation energy comparable to the energy of a single bond 
in the solid. 

The chemical potential p[R(z)] of a surface atom can 
be computed from the energy change due to the local ad- 
dition of the volume V a of an atom to the system. Within 
the Born-Oppenheimer approximation, and assuming the 
electrons act as an incompressible fluid [1, [12, 0], the 
chemical potential is given [3| by the functional deriva- 
tive p[R(z)} = V a /{2nR) ■ 5tt e /5R{z). Starting from Eqs. 
([1]) and J2]), one obtains 



fi[R(z)] = pa 



V a ( 2a s dC[R(z)} 8V t 



shell 



2irR 



V 1 + (d z Rf 



dR 



(5) 



LOsV a is the bulk chemical potential. Here 

Rd 2 R 



is the local mean curvature 



where po 

dC[R(z)]=ic(l- Tmw ^ 

of the wire, and arises from the functional derivative of 
the surface term in Eq. (p}. The chemical potential of a 




FIG. 1: Electron-shell potential V 3 heii(R) (top), chemical po- 
tential p C yi(R) (middle), and stability coefficient a(R) (bot- 
tom) for a cylindrical wire, as a function of its dimensionless 
radius kpR, where fef is the Fermi wavevector. (The radius 
range is limited to that corresponding to the simulations pre- 
sented here. For a graph over a more extended range, see Refs. 
[H, 13. fTH) Vertical dashed lines in the bottom panel mark the 
positions of conductance channel openings, which drive the 
instability [ll| . The top axis shows the conductance values of 
linearly stable cylinders in units of the conductance quantum, 
Go = 2e 2 /h. 
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cylinder /j, cy i(R) = u\R(z) = R] is plotted as a function 
of the radius R in the middle panel of Fig. [TJ 

The appropriate boundary conditions when simulating 
a cylinder, which in a real system is connected to larger 
leads [3, 9], have been shown to be Neumann boundary 
conditions, d z R — at both wire ends 0, 

Note that, despite the apparent simplicity of Eqs. ((3])- 
as written above, the resulting partial differential 
equation for R(z,t) is fourth-order in z-derivatives and 
highly nonlinear. Its classical coun terp art, with V s heii = 
0, has been extensively studied [25j, |26|, |27| . It was shown 
to have stationary states corresponding to shapes of con- 
stant mean curvature: The sphere, the cylinder and the 
unduloid of revolution, the latter being always unsta- 
ble. The addition of the electron-shell-potential term, 
of quantum-mechanical origin, is responsible for a rich 
dynamics, discussed extensively in Ref. [16j | and in this 
paper, as well as for the stabilization of the unduloid as 
a stationary state [3j, |9| . Without this term, any cylinder 
longer that its perimeter is unstable toward the Rayleigh 
instability and breaks apart into spheres. In that case, 
the system exhibits front propagation [28| , but there is no 
saturation of the instability, and the radius dependence 
of the dynamics is trivial. Furthermore, nonlinear effects 
do not seem to influence the front propagation. 



Summary: Linear stability analysis 

A linear stability analysis of cylindrical wires has been 
performed within the NFEM 0, EH H El . The change 
in the energy ((JJ) due to a radius perturbation AR — 
b q enp(iqz) was found to be, 

Att e = L^a(i? ;g)N 2 , (6) 

9 

so that the sign of a(Ro;q) determines the linear sta- 
bility of a wire of radius Rq towards a perturbation of 
wavevector q. Hence a(Ro; q) has been named the stabil- 
ity coefficient. It was further found that the global linear 
stability is essentially determined by the long wavelength 
limit [fJIIIOjI a(Ro) = a(Ro;q = 0), with 

, ( 2na s d 2 V she ii 1 dV she u\ 

a{Ra)= {-— + ^lR^-R^lR-) R ^ (?) 

Instabilities were found to result from a transverse 
eigenenergy of the wire crossing the Fermi energy Ep, 
thus closing or opening a conduction channel [3, 
These thresholds are marked by vertical dashed lines in 
the plot of a(R) in the bottom panel of Fig. [TJ Cylinders 
with radii in the vicinity of "magic" radii — corresponding 
to minima of the shell potential Vsheii i see Fig. [JJ — are 
linearly stable, while wires close to maxima of V s h e ii are 
unstable. 

A linearized dynamical theory [12| shows that unstable 
wires develop an exponentially growing instability with a 



well-defined wavelength A = 2%/q max , corresponding to 
the maximally unstable mode, such that a(i?o! Qmax) is 
extremal. This instability was argued to saturate, and 
eventually lead to a phase separation of the wire into 
thick and thin segments of stable radii [12j. Simulations 
using the full dynamics Eqs. ©"(O have confirmed this 
0], and shown that the phase separation occurs via a 
complex dynamics involving kink interactions and anni- 
hilation [l6|]. In this paper, I analyze the growth and 
propagation of such instabilities in greater detail. 



II. INSTABILITY PROPAGATION 

Results on the evolution of an unstable cylinder, in- 
cluding the growth of a single-wavelength perturbation 
followed by phase separation, were briefly presented in a 
Letter Q. In this section, I provide more details about 
simulations of the growth and propagation of perturbed 
linearly unstable wires as a function of the wire radius 
Rq. An initially localized perturbation is found to grow 
exponentially and propagate with a well-defined front, 
which is analyzed in detail. A theoretical analysis of the 
front propagation is provided in Sec. IIIIl 




FIG. 2: Instability propagation for a wire of radius /cf-Ro = 
11.65: The main panel shows <f>(z) = R(z) — Ro at equidistant 
times (shifted vertically for clarity), illustrating a typical time 
evolution of the instability propagation. In particular, the 
coherent pattern left behind the front can be observed. The 
r.h.s. axis gives the dimensionless evolution time r for the 
corresponding curves. The inset shows |0(z)| at a given time, 
r ~ 180, on a logarithmic scale, illustrating its exponential 
decay ahead of the front. 

The initial condition for these simulations is an un- 
stable cylinder of radius Ro, i.e. a radius such that 
a(R ) < 0. Two lengths, k F L = 300 and 400, have 
been considered. A Gaussian perturbation of amplitude 
UfSR = 0.01, localized at the left boundary, z = 0, 
is added to trigger the instability. Results have been 
checked to be independent of the type, amplitude, and 
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extension of the initial perturbation, as long as it is suf- 
ficiently localized 

After a brief initial incubation period, during which 
the amplitude of the perturbation decreases, it grows 
exponentially and propagates along the wire. A well- 
defined front forms and moves at a constant velocity, 
as illustrated in Fig. 03 where the wire perturbation 
4>(z, t) — R{z,t) — Rq for a wire with kpRo — 11.65, 
is plotted at equidistant times r, with a vertical shift 
proportional to r. Behind the front, a coherent pattern 
consisting of a single- wavelength radius oscillation forms 
and saturates at an amplitude kp/S.(f) ~ 0.1—1, whose 
value depends on the details of the shell potential V s heii 
(cf. Fig. [T]) in the vicinity of the unstable radius Rq con- 
sidered. 

A logarithmic plot of (the absolute value of) the wire 
perturbation |$(z,t)| (see inset of Fig. [3]) clearly shows 
its exponential decay ahead of the front, whose position 
Zf is defined as the point where the perturbation ampli- 
tude reaches a certain threshold value c [29J. Tracking 
Zf as a function of time, the front velocity v can be ex- 
tracted, and is found to be constant after an initial period 
when it is influenced by boundary and initial conditions. 
Similarly, the wavevector q and decay length £ of the per- 
turbation ahead of the front can be extracted. Both are 
found to be constant after an initial decrease during early 
evolution. 

In the simpler cases, further evolution of the wire be- 
hind the front occurs on a timescale orders of magni- 
tude longer than the front propagation. It involves an- 
nihilation of kink/antikink pairs — upward and downward 
steps in R(z), whose dynamics has been described in Ref. 
[16| — and eventually leads to phase separation into thick 
and thin segments of stable cylinders [H, [HJ • 

In other cases however, nonlinear effects quickly take 
over the dynamics behind the front, and start pushing it. 
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FIG. 3: log \tj>(z, t)/0o| at equidistant times r, shifted verti- 
cally for clarity, for a wire of radius UfR = 12.28, showing 
an invasive mode. All curves are normalized by (f>o = 4>( z = 
0, r — tq), where To is the time of the bottom curve. 



This modifies its velocity, as well as the decay rate, and 
possibly the wavelength ahead of the front. An example 
of such behavior is shown in Fig. 02 where log \<f){z, r)| 
is plotted at equidistant times (each curves shifted ver- 
tically) for a wire with radius kpR = 12.28. One can 
clearly see a second front with a shorter decay length 
that progressively takes over the initial front. Such a 
front will be referred to as a nonlinear front, while the 
initial, slower front will be called linear. 

The velocity, decay length and wavevector of both 
fronts can be extracted from the simulations. The three 
quantities are plotted as open circles in Fig. 0] for the lin- 
ear front, while filled light (green) circles are used for the 
nonlinear front, when it is detected. All three quantities 
are found to depend non-trivially on the wire radius Rq. 
They reach local extrema simultaneously for values of Rq 
that correspond to minima of the stability coefficient a 
(Fig. CQand Eq. 0), i.e., maximally unstable wires. The 
front velocity v and wavevector q vanish, while the decay 
length £ diverges, at the stability boundary R c , where 
a(R c ) = 0. Roughly, it seems that a nonlinear front 
appears for wires with a(Ro) > Q.lEpkp. 



III. FRONT PROPAGATION ANALYSIS 

In this section, I provide a more detailed analysis of the 
instability propagation, showing in particular that the 
evolution of the linear fronts — marked by open circles 
in Fig. [4] — derives from the linearized dynamics, hence 
the name, and is a typical example of front propagation 
into an unstable state. Following ideas from the linear 
marginal stability analysis of Refs. [l|, |3fJ , much of the 
front dynamics can be understood from the linearized 
evolution equation for the perturbation <fi(z,t), which is 



2tt.Ro ■ 



V a dt 



27R7 s i? — -a(R )-^ 



dz 2 



(8) 



Assuming a front of the form <j){z,t) = cxp[i(a;i — kz)\, 
with us and k complex quantities, its dispersion relation 



w{k) 



E, 



k z \k' + 



a(Ro) 

27TCT s i? 



(9) 



is derived. It has been argued [30] that the front wavevec- 
tor k* is such that the front is marginally stable, and thus 
satisfies the conditions [l[ 



dto 
Im dk 



0. 



Imcj(fc*) 
Im k* 



dk 



(10) 



This yields the front wavevector k* = q — 




(11) 



where the real part q determines the wavevector of the 
pattern left behind the front, and the inverse £ of the 
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FIG. 4: (color online) Instability wavevector q (top), decay 
length £ (middle) , and front velocity v (bottom) as a function 
of unstable wire radius Ro- The solid lines give the results 
of the linear theory, Eqs. and (|13|l . and the circles are 
results extracted from the full dynamical simulations (Error 
bars are smaller than the symbols, and therefore not plotted.) 
Open circles are used for the linear front, while filled (green) 
circles correspond to nonlinear fronts. The velocity in the 
bottom panel is given in units of vo = Xfu>o- 



imaginary part corresponds to the decay length ahead 
of the front. They are respectively plotted as solid lines 
in Fig. [4] (top two panels), and compare well with data 
extracted from the full nonlinear dynamics (open circles) 
of initial fronts. Small deviations are observed for large 
values of £, but are a result of the difficulty of determining 
a large decay length in a finite system. 
It is clear from Eq. (fTTj) that the product 



'3(3 + 77) 
V7-1 



3.2, 



(12) 



is universal and does not depend on any wire parame- 
ters. This has been verified within numerical accuracy 
for all wire radii, as well as for two different values of a s 
(corresponding to Au and Na). 

The front velocity v — Im[w(fc*)]/Im k* is readily ex- 
tracted from Eq. (fTTj) . and is found to be 



uj a o s V 2 a / 17 + 7^7 / a(R ) 



6E F 



1T(T°R 



3/2 



(13) 



This result is plotted in Fig. [4] (bottom panel) as a solid 
line, and compared with linear-front velocities obtained 



from the dynamical simulations (open circles). Agree- 
ment between the linear theory and full nonlinear dy- 
namics is very good, showing that the front propagation 
is indeed governed by the linear dynamics. 
Combining Eqs. (fTTj) and (fl3jl. one gets 



(14) 



v oc a s q 



so that the front velocity can be determined from the 
pattern wavelength. This relation also holds within nu- 
merical accuracy in simulations for Au and Na. 

Finally, as a(R) oc ±(i? — R c ) at the stability bound- 
aries, where a(R c ) = 0, Eqs. (fTTj) and (flU)) provide the 
scaling of v, q, and £ as \R — R c \ — > 0, which are 



\R - Rc 



3/2 



\R-R, 



1 1/2 



and 



£ ~ \R - R, 



1-1/2 



(15) 



(16) 



(17) 



In cases where nonlinear effects take over the dynamics 
behind the front, the marginally stable front satisfying 
Eqs. (fmfTTj) forms and propagates for a while, but at 
some point it is "invaded" by a faster front. This "in- 
vading" front has a larger q and a smaller £, and satisfies 
neither Eq. (JT3J) , nor Eq. (fT4"|) . This is consistant with the 
nonlinear marginal-stability mechanism discussed by van 
Saarloos [3l| where, for some range of parameters, the 
front becomes unstable to an "invasion mode." In this 
case, the evolution depends on the full nonlinear dynam- 
ics, Eqs. J5HS1), and the front is "pushed" by the invading 
mode, rather than "pulled" by the exponentially growing 
instability [1]. The new front speed, wavevector and de- 
cay length seem consistent with the analysis of Ref. [3l|. 
Although the cause of the instability of the front to the 
invasive mode is not clear, it seems to be related to the 
existence of a "quasi-stable" wire, i.e. a wire for which 
a(R), though still negative, is relatively small. 



IV. DISCUSSION AND CONCLUSIONS 

Experimental verification of the dynamics of an unsta- 
ble wire may be, at least in part, possible. An unstable 
cylinder can be prepared using a potential bias: As the 
magic radius intervals vary with a bias applied along the 
wire j321 | . a stable cylindrical wire can be prepared at 
a finite bias V, and its stability modified by suddenly 
switching the bias off. If the bias and wire radius are 
chosen appropriately, the wire becomes unstable, and its 
dynamics can be observed using, for example, transmis- 
sion electron microscopy Q. An alternative way of ob- 
taining the same result is to stretch the wire abruptly, 
so that it deforms elastically into an unstable cylinder. 
In both cases, the connection of the wire to macroscopic 
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leads acts as a localized perturbation at the end of the 
wire triggering the instability propagation. 

Actual observation of the front propagation would, 
however, be difficult as wire imperfections are likely to 
trigger the instability at several places along the wire si- 
multaneously. In addition, the experiment needs to be 
conducted at room temperature for the metal to be soft 
enough to allow surface diffusion, so that thermal fluc- 
tuations are likely to have the same effect. Furthermore, 
the predicted radius oscillations are of the same size as 
the atomic granularity, and the two may thus be difficult 
to distinguish. 

The evolution of unstable wires, discussed in Sec. [TTl 
provides an interesting system where ideas developed in 
the context of front propagation into an unstable state 
0, [3(| HH can be successfully tested, as shown in Sec. 
IIIII The instability is found to grow exponentially ahead 
of a well-defined front, leaving a coherent pattern behind 
it. In most cases, the front is found to be pulled by the 
instability growth ahead of it. Its propagation is thus 
governed by the linearized evolution equation ([5]). Sim- 



ple expressions for the front velocity v and decay length 
£, as well as the pattern wavevector q have been derived, 
and found to be in good agreement with numerical sim- 
ulations using the full non- linear dynamics. 

In other cases, the linear evolution equation © fails to 
explain the full front dynamics. The pulled front is found 
to be invaded by a faster instability that pushes the front 
at its higher velocity. In this case, the front parameters 
depend on the full non-linear dynamics. The presence 
of an invading mode seems to correspond to quasi-stable 
wires, for which a(i?o) is still negative, but small. 
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